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ABSTRACT: The quenched approximation for QCD is, at present and in the foreseeable 
future, unavoidable in lattice calculations with realistic choices of the lattice spacing, vol- 
ume and quark masses. In this talk, I review an analytic study of the effects of quenching 
based on chiral perturbation theory. Quenched chiral perturbation theory leads to quan- 
titative insight on the difference between quenched and unquenched QCD, and reveals 
clearly some of the diseases which are expected to plague quenched QCD. 
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1. Introduction 



The lattice formulation of QCD has proven to be a powerful tool for computing QCD 
quantities of direct phenomenological interest, such as hadron masses, decay constants, 
weak matrix elements, the strong coupling constant, etc. (For reviews see for instance refs. 
[1,2], or the proceedings of Lattice 93 [3].) 

In order to perform such computations numerically, one obviously needs to consider a 
system with a finite number of degrees of freedom, which is accomplished by putting 
lattice QCD in a finite box. This box is then hopefully large enough to accomodate the 
physics one is interested in without serious finite size effects. This leads to the requirement 
that the Compton wavelength of the particles of interest is sufficiently smaller than the 
linear dimension of the box, i.e. the mass has to be large enough for the particles to fit in 
the box. 

In order to have a small enough lattice spacing, small enough masses (in particular for 
the pion) and a large enough box size, one needs a large number of degrees of freedom in 
a numerical computation. It turns out that for QCD with realistic choices of the lattice 
spacing a, volume V and the quark masses (in particular the light quark masses), the 
presently available computational power is not adequate. The most severe problem comes 
from the fermion determinant, the logarithm of which is a very nonlocal part of the gluon 
effective action (specially for light quark masses). This nonlocality slows down the Monte 
Carlo algorithms dramatically. 

In order to circumvent this problem, most numerical computations in lattice QCD have 
been done in the quenched approximation, in which one simply sets the fermion determi- 
nant equal to one [4]. This amounts to ignoring all fermion loops which occur in QCD 
correlation functions (except those put in by hand through the choice of operators on 
the external lines). While some handwaving arguments exist as to why this might not 
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be unreasonable, the quenched approximation does introduce an uncontrolled systematic 
error. Since the effect of a fermion loop is roughly inversely proportional to its mass, this 
error is expected to be particularly severe for quantities involving light quarks. It appears 
therefore that chiral perturbation theory ChPT maybe a useful tool for investigating the 
difference between quenched and unquenched ("full") QCD. 

In this talk, I will review a systematic approach to the study of the quenched approximation 
through ChPT [5,6,7,8]. There are two reasons why ChPT is useful in this context: 

• It turns out that ChPT can be systematically adapted to describe the low energy sector 
of quenched QCD [6]. It will therefore give us nontrivial, quantitative information on the 
difference between quenched and full QCD. 

• ChPT describes the approach to the chiral limit, and can be used for extrapolation of 
numerical results to small masses and large volumes. If these results come from quenched 
computations, one will of course need a quenched version of ChPT. (For finite volume 
ChPT, see refs. [9]. For quenched finite volume results, see refs. [6,7].) 

In this review, I will concentrate on the first point. I will first show how ChPT is developed 
for the quenched approximation, and then use it for a quantitative comparison between 
full and quenched QCD. The quantities that I will discuss are fx/fn [6,8] and the octet 
baryon masses [10]. 

I will then address a number of theoretical problems that arise as a consequence of quench- 
ing. That such problems arise is no surprise, as quenching QCD mutilates the theory quite 
severely. It is however quite instructive to see what the detailed consequences are. 

2. Systematic ChPT for quenched QCD 



In this section I will outline the construction of a chiral effective action for the Goldstone 
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boson sector of quenched QCD [6]. I will first introduce the formalism, and then show how 
it works in some examples. For early ideas on quenched ChPT, see ref. [11,5]. 

We will start from a lagrangian definition of euclidean quenched QCD. (We will restrict 
ourselves entirely to the euclidean theory which can be defined by a pathintegral. Hamil- 
tonian quenched QCD presumably does not exist.) To the usual QCD lagrangian with 
three flavors of quarks q a , a = u,d, s, we add three ghost quarks q a with exactly the same 
quantum numbers and masses m a , but with opposite, bosonic, statistics [11]: 

^quarks = Yl (P + ™a)?a + X^a(^ + m a )q a , (1) 
a a 

where jp is the covariant derivative coupling the quark and ghost quarks to the gluon field. 
The gluon effective action produced by integrating over the quark- and ghost quarkfields 
vanishes, since the fermion determinant of the quark sector is exactly cancelled by that of 
the ghost sector. Note that the ghost quarks violate the spin-statistics theorem. 

We will now assume that mesons are formed as (ghost) quark - (ghost) antiquark pairs 
just like in ordinary QCD. This is basically equivalent to the notion that it is the dynamics 
of the gluons which leads to confinement and chiral symmetry breaking. The Goldstone 
particle spectrum of quenched QCD will then contain not only qq, but also qq, qq and qq 
bound states. We will denote this 36-plet by 

<$>=(<!> x!) ^(qq qq) (2) 

\x <t> J \qq qq) 

Note that the fields x an d X 1 describe Goldstone fermions. 

The quenched QCD lagrangian eq. (1) for vanishing quark masses has a much larger 
symmetry group than the usual U{3)i x U(3)r flavor group; it is invariant under the 
graded group U(3\3)l x U(3\3)r [6], where U(3\3) is a graded version of U(6) since it 
mixes the fermion and boson fields q and q. Writing an element U of U(S\S) in block form 

as 

V=(* ° B ), (3) 
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the 3x3 matrices A and B consist of commuting numbers, while the 3x3 matrices C and 
D consist of anticommuting numbers. 

We can now construct a low energy effective action for the Goldstone modes along the 
usual lines. We introduce the unitary field 

E = exp(2i$/f), (4) 

which transforms as E — > U^HU^ with and Ur elements of U(3\3). Because we are 
dealing here with a graded group, in order to build invariants, we need to use the supertrace 
str and the superdeterminant sdet instead of the normal trace and determinant, with [12] 

str(U) = tr(A) - tr(B), 
sdet(U) = exp(str log (U)) = det(A — CB~ 1 D)/det(B). (5) 

To lowest order in the derivative expansion, and to lowest order in the quark masses, the 
chiral effective lagrangian consistent with our graded symmetry group is 

f 2 

£o = J —str{d^d^) - v str{MT, + MTt), (6) 
8 

where M. is the quark mass matrix 



/ and v are bare coupling constants which are not yet determined at this stage. 

The symmetry group is broken by the anomaly to the smaller group 
[SU(S\S)l x SU f>](s)U (1) (the semidirect product arises as a consequence of the graded 
nature of the groups involved; the details are irrelevant for this talk). SU(3\3) consists 
of all elements U G f7(3|3) with sdet(U) = 1. The anomalous field is <E>o = (7/ — fj')/y/2, 
where the relative minus sign comes from the fact that in order to get a nonvanishing tri- 
angle diagram, one needs to choose opposite explicit signs for the quark and ghost quark 
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loops, due to the different statistics of these fields, rj is the field describing the normal rf 
particle, while fj' is the ghost rj' consisting of ghost quarks and ghost antiquarks. We will 
call the field $o the super-?/ field. The field $o ^ log ^ = log ^ is invariant under 
the smaller symmetry group, and we should include arbitrary functions of this field in our 
effective lagrangian. Following ref. [13], the correct chiral effective lagrangian is 

C =yi($ >M<V^ Et ) - V 2 {$ )str{MX + M£ f ) 

+ y ($o)+V5($o)(9 M $o) 2 , (8) 

where the function multiplying i str{MY> — MT^) can be chosen equal to zero after a field 
redefinition. This lagrangian describes quenched ChPT systematically, as we will show 
now with a few examples. 

For our first example, let us isolate just the quadratic terms for the fields rf and fj' , choosing 
degenerate quark masses for simplicity. We expand 

Vi($ ) = y + ---> 

y 2 ($ ) =u + ..., 

V ($ ) = constant + /J 2 $o + --- 5 ( 9 ) 
l/ 5 ($ ) =« + ..., 

and obtain 

=^(<V/) 2 - ^(<V/') 2 + \a(d^> - d^fj'f 

+ ^m 2 (r/) 2 - \m\{fif + \f(r{ ~ v'f + • • • , (10) 

where m% = 8mv/ f 2 . The relative minus signs between the rj' and fj' terms in eq. (10) 
come from the supertraces in eq. (8), and are related to the graded nature of the chiral 
symmetry group of quenched QCD. 
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The inverse propagator in momentum space, 

(P 2 +™1)(1 _° 1 )+(^ 2 + ap 2 )(_ 1 1 (11) 

clearly cannot be diagonalized in a p independent way, which is quite different from what 
one would expect from a normal field theory! Treating the y? + up* term as a twopoint 
vertex, one can easily show that the repetition of this vertex on one meson line vanishes, due 
to the fact that the propagator matrix ^ J ~~ 1 ) multiplied on both sides by the vertex 

matrix ^ ^ ^ gives zero. This result coincides exactly with what one would expect 

from the quark flow picture for 7/ propagation, as depicted in fig. 1. The straight-through 
and double hairpin contributions do not contain any virtual quark loops, and are therefore 
present in the quenched approximation. All other contributions should vanish because 
they do contain virtual quark loops, and this is exactly what happens as a consequence 
of the (admittedly strange) Feynman rules for the propagator in the rj'-fj' sector! This 
propagator is given by the inverse of eq. (11) and reads 



1 \ f 2 + ap 2 ( 1 1 



p 2 + ml V° - 1 ) (;p 2 + m2)2 \l 1J> 
in which the two terms correspond to the two first diagrams in fig. 1. 



(12) 



Figure 1. Quark flow diagrams for the 7/ propagator in full QCD. 

From eq. (12) we learn several things. First, because y? , which in full ChPT would 
correspond to the singlet part of the rj mass, appears in the numerator, we need to keep 
the 7/ (and its ghost partner) in quenched ChPT: it cannot be decoupled by taking y? 
large. Second, this "propagator" is definitely sick, due to the double pole term. It should be 
stressed here that this double pole term is an unescapable consequence of quenched QCD, 
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and does not result from our way of setting up chiral perturbation theory. In the case of 
nondegenerate quark masses, this double pole also shows up in the 7r° and rj propagators, 
due to mixing with the rf . I will return to these strange properties of quenched QCD in 
section 4. 



As a second example, we will consider the (charged) pion selfenergy at one loop, again with 
degenerate quark masses. I will set a = for simplicity. At one loop, the pion selfenergy 
only contains tadpoles, with either or x lines (cf. eq. (2)) on the loop. Also, on the 4> 
loop, one can have an arbitrary number of insertions of the vertex (j 2 if the internal <p line 
is an SU(3) singlet. These various contributions are drawn in fig. 2, where a solid line 
denotes a <fi line, a dashed line denotes a x nne ? an d a cross denotes a y? vertex. One finds 
that the diagrams with the 4> and x lines on the loop without any crosses cancel, and then, 
of course, that the diagrams with more than one cross vanish, using our earlier result for 
the ij'-fj' propagator. We are left with only one term, and the result is 



The pion selfenergy is logarithmically divergent, but the origin of this divergence is com- 
pletely different from those that arise in the unquenched theory, as it is proportional to 
fx 2 . One can easily convince oneself that the diagrams in fig. 2 which cancel or vanish 
correspond to diagrams with virtual quark or ghost quark loops in the quark flow picture. 
(For early discussions of the quenched pion selfenergy in the quark flow picture, see refs. 



Before I go on to look at some quantitative results, I would like to discuss one aspect of 
the chiral expansion in quenched ChPT. The chiral expansion is basically an expansion 



Figure 2. One loop pion selfenergy in quenched ChPT. 




(13) 



[11,5].) 
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in the pion mass (see e.g. ref. [14]). However, as we have argued above, in quenched 
ChPT there is unavoidably another mass scale, namely the singlet part of the rf mass, 
y?. For our expansion to be systematic as an expansion in the pion mass, we would 
have to sum up all orders in y? , at a fixed order in the pion mass. This is clearly a 
formidable task. In order to avoid this complication in a systematic way, we can think of 
// 2 /3 (which turns out to be the natural parameter as it appears in the chiral expansion) 
as an independent small parameter. To check whether this makes any sense, one may note 
that from the experimental value of the 7/ mass one obtains a value /U 2 /3 ~ (500 MeV) 2 , 
which is roughly equal to the kaon mass squared, m? K - Of course, for quenched QCD 
the parameter y? need not have the same value, after all quenched QCD is a different 
theory. A lattice computation of this parameter [15] gives ^quenched/^full ~ 0-75. (a can 
be estimated from rj-rj' mixing, and is very small.) Finally, one may also note that both 
y? and a are of order 1/N C , where N c is the number of colors [16]. I will return to this 
point in section 4. 

3. Quantitative comparison of quenched and full ChPT 

Let us first consider the quenched result for the ratio of the kaon and pion decay constants 
fx and f n [6,8]. I will set a = and take m u = = m: 

+ (m s -m)L. (14) 

L is a certain combination of "low energy constants" [13]. Since this constant is a bare 
parameter of the quenched chiral lagrangian, the result eq. (14) is not directly comparable 
to the equivalent result for the full theory. In other words, in order to compare quenched 
and full QCD, we have to consider quantities which are independent of the bare parameters 
of the effective action. (Alternatively, we would need to extract the values of the bare 
parameters from some independent measurement or lattice computation, in this case, we 
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would need independent determinations of L in both quenched and full QCD.) In the full 

1/3 4/3 

theory, (frjfn )/ fx * s such a quantity [13], but in the quenched theory this quantity is 
not well defined, due to the double poles which occur in the propagators of neutral mesons. 

We will therefore choose to consider a slightly different theory, in which sufficiently many 
decay constants referring only to charged (i.e. off-diagonal) mesons are present [8]. This 
theory is a theory with two light quarks m u = = m and two heavy quarks m s = m s i = 
m' . This theory contains a ud pion n, an s's pion n' and a us kaon K, with mass relation 



2 1/2,2 
m K = 2^ + m ^ 



(15) 



One can show that the ratio fxl\/ firfn' ^ s independent of the low energy constants L. 
For the quenched theory we find 



fx 



\ 1— loop 



V 1 fitf-K 1 J quenched 

whereas in the full theory 

1— loop 



1 + 
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ml + m 2 ^ 
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2(7774,-7774) \m 



- 1 



(16) 
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K 
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(17) 



Note again that the logarithms in the quenched and unquenched expressions are completely 
different in origin. 

We may now compare these two expressions using "real world" data, where we'll determine 
the value of the it' mass from the mass relation eq. (15). With m^ = 140 MeV, m^ = 
494 MeV and ^ 2 /3 = 0.75 x (500 MeV) 2 we find 



f K 



fK 



1— loop 

quenched 
1— loop 



1.049, 



1.023, 



(18) 



. \/J~f-' / ( ull 

a difference of 3%. If we choose ^2/3 = (500 MeV) 2 , we find a difference of about 4%. 
This difference is small. Note however, that this is due to the fact that for this particular 
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ratio, ChPT seems to work very well, both for the full and the quenched theories. If one 
only considers the size of the one loop corrections (the numbers behind the decimal point), 
the quenched and full results are very different. It is also possible, and in fact not unlikely, 
that part of the difference between the full and quenched theory gets "washed out" by the 
fact that we are considering a "ratio of ratios" . It follows that the relative difference is a 
lower bound on the difference between the quenched and full values of the decay constants. 
For another quantity for which the difference between quenched and full ChPT has been 
calculated, see ref. [8]. 

Next, I will review some very recent work on baryons in quenched ChPT by Labrenz and 
Sharpe [10]. They calculated the one loop corrections to the octet baryon mass coming 
from the cloud of Goldstone mesons. They employed an effective lagrangian for quenched 
heavy baryon ChPT, constructed using the same techniques as described in section 2. In 
the case of degenerate quark masses, the result for the nucleon mass is 



In this equation, m, D, F, bp,, bp and 7 are bare parameters which occur in the baryon 
effective action, m is the bare "average" octet mass, D and F are the well known baryon- 
meson couplings, bp, and bp are low energy constant which arise as a consequence of 
renormalization (see for instance refs. [17,18]). 7 is a new coupling which occurs because 
of the unavoidable presence of the super-?/ in the quenched approximation. The term 
proportional to /j, comes from a diagram with a cross on the <f> internal line, i.e. an insertion 
of the p? twopoint vertex. Note that in this case there are also one loop corrections not 
involving \j? which survive the quenched approximation, in contrast to the pion selfenergy, 
eq. (13), or fx/ fn, eq. (14). The authors of ref. [10] then calculated the coefficients using 
full QCD values for the various parameters (from ref. [18]). With a = and 7 = (7 = 



m N =m-—(D- ZFff^mv + 2(b D - 36 F )mJ 




7T 



+ ^(D-3F)(2D + 3 1 ) + ^a(D-3F) 2 




(19) 
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is consistent with available information, which however is limited [19]), they obtained 

m N = 0.97-0.5— m n + SAml- 1.5mJ, (20) 
with 5 = fi 2 /(24ii 2 f%) and 5 w 0.2 for the full theory. 

In ref. [10], eq. (19) was also compared to recent numerical results from ref. [20]. These 
data are presented in fig. 3, where the scale a~ l = 1.63 GeV is set by f n [10]. If one 
calculates the coefficients in eq. (19) by "fitting" the four data points, one finds 

m N = 0.96 - I.Otott + 3.6ml ~ 2.0mJ. (21) 

This is only four data points for four parameters, and the "fit" is quite sensitive to for 
instance an additional vn\ term. From the agreement between eq. (20) and eq. (21) 
it appears that it is reasonable to apply ChPT to the results of ref. [20]. Note that 
the individual terms in eq. (20) are quite large for the two higher pion masses in fig. 
3 (this is not unlike the case of unquenched ChPT). From fig. 3 it is also clear that 
(m/ /7r)quenched 7^ (^//nOfull because of the term linear in eq. (19), which is absent in full 
ChPT. 

Labrenz and Sharpe then went on to consider octet mass splittings. In order to remove 
effects which can be accomodated by a change of scale, they calculated the ratios 

Tfl ' 

Rij = —, U = iV,A,£,£ (22) 
mj 

in quenched ChPT, and compared these with similar ratios obtained from ref. [18]. They 
assumed that all bare parameters in the equations for the octet masses (for explicit ex- 
pressions, see their paper) are equal in the full and quenched theory, and then calculated 
the ratios 

D quenched 

T ij = ofull • ( 23 ) 
R ij 
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With the assumption that the bare parameters of the quenched and full theories are equal, 
bp and bp drop out of the ratios, and with 7 = 0, a = and D and F equal to their full 
QCD values, they obtain 

r^ N = 1 + 0.19(5/0.2) + 0.13 = 1.31[1.27] for 5 = 0.2[0.15], 

r~N = 1 - 0.46(5/0.2) + 0.43 = 0.97[1.09] for S = 0.2[0.15], (24) 

r AN = 1 - 0.39(5/0.2) + 0.26 = 0.87[0.97] for 5 = 0.2[0.15]. 

(The choice 5 = 0.15 corresponds to the value reported in ref. [15].) 
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Figure 3. The nucleon mass from the lattice [20] (copied from ref. [10]). The curve is 
from a fit to the form rajy =rn + am^ + bm^ + cm\ . 

From this, one would conclude that one can expect errors from quenching of at least 20% 
in the octet splittings. These differences between the quenched and full theories cannot be 
compensated for by a change in scale between quenched and full QCD. 

At this point I would like to comment on the above mentioned assumption that was used 
in order to obtain eq. (24). Let us consider in particular the parameters bp and bp. They 
correspond to higher derivative terms in the baryon-meson effective action, and are needed 
in order to absorb the UV divergences which arise at one loop in ChPT. Since the size 
of these divergences is in principle different between the full and quenched theories, one 
expects that ^quenched and ^quenched can be different from b D j u n and ftjrfuu- If we 
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want to proceed without assuming that the quenched and full 6's are equal, we have to 
consider ratios of quantities independent of the parameters bp, and bp. The situation is 
essentially the same as in the case of fx/ fir- From the available results [10], only one ratio 
independent of bp and bp can be formed: 

X = ^4-. (25) 

If we expand X in the Goldstone meson masses using ChPT, X — 1 measures the deviation 
from the Gell-Mann-Okubo formula (cf. ref. [18] for the full theory). 

Setting rriTY = keeping only mx as in ref. [10], one finds 

3 

^quenched = 1 " 1.1046^^2 ~ 3i?2 ) + L33335 7f= ~ 3i?2 ) ' 

X m = 1 - 0.4125^^2 (d 2 - 3F 2 ) . (26) 

(The parameter 7 drops out of this particular combination and we again take a = 0.) 
These quantities still depend on the other bare parameters, D, F and m. Again, they 
could be different in the quenched and full theories, and I will leave the quenched values as 
free parameters. Substituting m K = 495 MeV, f n = 132 MeV, D in y[ = 0.75, F full = 0.5 
and mf u n = 1 GeV [18] finally gives 



^=^ = 1-0.0214 + 
^full 



5 

0.293 0.306 

0.2 



/ r)2 3 

V / quenched 

^quenched/ 1 GeV 



(27) 



For any reasonable values of m and 6, and for (D 2 — 3-F 2 ) qucnclied not too far from its full 
theory value of —0.1875, the difference between the quenched and full theories as measured 
by the ratio ^ qu enched/^full i s n °t more than a few percent. Of course the same comment 
that applied in the case of fx/ fir applies here, that part of the difference may have been 
washed out by taking "ratios of ratios" . Summarizing, the conclusion of this analysis seems 
to be that the error from quenching for octet baryon masses is at least a few percent, and 
could be as much as 20%. 
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4. A sickness of quenched QCD 



Let us again consider the quenched result for fx/ fir, ec l- (14), as a function of the quark 
masses (using treelevel relations between meson masses and quark masses), 



M 1_loop _ 1 , kV3 

/ quenched 167T 2 /| 



m u + m s m s 

log 1 



+ (L — terms). 



_2(m s - m u ) m u 

From this expression it is clear that we cannot take m u — > keeping m s fixed, or, to 
put it differently, that if we take both m u and m s to zero keeping the ratio fixed, the 
limit depends on this ratio, and is not equal to one! This is quite unlike the case of full 
ChPT, where one can take any quark mass to zero uniformly, and deviations from SU (3) 
symmetry due to this quark mass vanish in this limit. Technically, the reason for this 
strange behavior is that there is another mass /x, which, as we argued before, cannot be 
avoided in quenched ChPT. This mass is related to the singlet part of the rf mass, and is 
not a free parameter of the theory. Even if we do not consider any Green's functions with 
rf external lines, this mass shows up through the double pole term in eq. (12) on internal 
lines. Because of the double pole, such contributions can lead to new infrared divergences 
in the m n — > limit. This problem with the chiral limit of quenched ChPT shows up in 
other quantities, such as meson masses and (ipip) [6,7,21]. 

A question one might ask is whether this problem is an artifact of one loop quenched 
ChPT [8]. For instance, if we would sum all contributions to the rf propagator, maybe the 
double pole term would become softer in the p — > limit, improving the infrared behavior 
of diagrams in which the double pole terms appear. Let us address this question in the 
chiral limit, m a = 0, where the problem is most severe. In the full theory, we can write 
the fully dressed rf propagator as 

Z{V) (28) 



p 2 + E(p) ' 

and define fj.'p(p) = H(p), which onshell is the rf mass in the chiral limit. Likewise, in the 
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quenched theory we can write the dressed r\ , r\ propagator as 




which defines Hq(p). To leading order in 1/N C , these two definitions of fi (p) should lead 
to the same result: 

We also believe that fJ/p{p = 0) is not equal to zero, since we expect the 7/ to remain a 
well-behaved, massive particle in the chiral limit. This implies, in sofar as we can rely on 
the large N c expansion, that A*q(0) 7^ 0, and that the double pole in eq. (12) is a true 
feature of the theory. 

While this argument is not very rigorous, I believe that the foregoing discussion implies 
that the chiral limit of quenched QCD really does not exist. This believe is futhermore 
supported by the following remarks: 

• Sharpe [7] has summed a class of diagrams in the case of degenerate quark masses for a 
very simple quantity (the pion mass), and found a result that is actually more divergent 
than the one loop result. 

• With nondegenerate quark masses there are many more diagrams that are infrared 
divergent in the chiral limit, and it is even less probable that resummation will improve 
the situation. 

• Any mechanism improving the infrared behavior would have to work for each divergent 
quantity. One expects that such a mechanism would be related to the double pole term in 
the rj' propagator, which created the problem in the first place. But this seems unlikely in 
view of the arguments given above. 

• The bare quark mass parameter appearing in the chiral effective action is not the same 
as that appearing in the (unrenormalized) QCD lagrangian. But one can argue that the 
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two bare quark masses should be analytically related, and the infrared problem is not just 
a problem of quenched ChPT, but of quenched QCD. 

5. Conclusion 

The quenched approximation leads to an unknown systematic error in all lattice calcula- 
tions that employ this approximation. It would of course be very nice to have a parameter 
that interpolates between full and quenched QCD, and in principle the quark masses could 
play such a role, since one expects that quenched QCD corresponds to full QCD with 
very heavy quarks. One would have to distinguish here between valence and sea quark 
masses, and it is the sea quark mass that would play the role of such a parameter. This 
distinction can indeed be made by considering so-called partially quenched theories [22], 
but no practical scheme to implement this idea is known. 

Quenched QCD can be defined from a euclidean pathintegral as rigorously as full QCD. 
In this talk I have explained that euclidean quenched ChPT can be used as a tool for a 
systematic investigation of quenched QCD. Quenched ChPT does not quite accomplish a 
task equivalent to that of an interpolating parameter. Since the bare parameters appearing 
in the quenched and full chiral effective actions are not the same (as explained in section 3) 
one cannot directly compare quantities calculated in full and in quenched ChPT. However, 
one can calculate combinations of physical quantities which do not depend on the bare 
parameters, and in that case a direct comparison between quenched and full QCD is 
possible, as we demonstrated with an example involving meson decay constants. This 
makes it possible to estimate lower bounds on the differences which come from quenching; 
these estimates are dependent on the values of the meson masses, which can be taken to 
be the (known) independent parameters of the theory. For realistic values of these masses, 
such differences turn out to be of the order of a few percent for ratios of decay constants 
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and for baryon octet splittings. 

The disadvantage of this more conservative approach is that part of the difference maybe 
hidden, because these specific combinations of physical quantities maybe less sensitive to 
the effects of quenching than other quantities of interest. This is particularly clear in 
the example of baryon octet masses. In this case, a comparison based on the assumption 
that the bare parameters of the full and quenched effective theories are the same, lead 
to differences of up to 20% and more. Of course, it is not known to what extend this 
assumption is valid. 

The differences between the quenched and full theories become markedly larger for decreas- 
ing quark masses. This is due to the fact that new infrared divergences occur in quenched 
QCD, which do not have a counterpart in full QCD. These divergences lead to the nonex- 
istence of the chiral limit for quenched ChPT (as discussed in section 4). The origin of this 
phenomenon can be traced to the special role of the rf in the quenched approximation. 
In the quenched approximation, the rf is a Goldstone boson (it develops massless poles in 
the chiral limit), but an additional double pole term arises in its propagator, rendering it 
a "sick" particle. For nondegenerate quark masses this problem is also inherited by the n° 
and the r\. In section 4 I argued that the nonexistence of the chiral limit is a fundamental 
feature of quenched QCD. 

In principle therefore, the chiral expansion breaks down for quenched QCD. For very small 
quark masses, at fixed values of the singlet part of the rj' mass fj 2 , the expansion becomes 
unreliable. In order to make progress, one may take the expansion to be an expansion in 
// 2 /3 (which was shown to be roughly equal to m K phenomeno logically), with coefficients 
which are functions of the quark mass. These functions then can be expanded in terms of 
the quark masses, sometimes leading to divergent behavior of the leading term {e.g. the 
one loop correction to fx/ fn)- If such divergences occur, the expansion is only valid for a 
range of quark masses which are neither too small, nor too large. It would be interesting 
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to see whether this point of view can be made solid. 

It is in principle interesting to study any quantity which is being computed in quenched 
lattice QCD in ChPT, for those quantities for which ChPT is applicable (meson masses, 
decay constants, condensates and the kaon B parameter have been calculated [6,5,8]). As 
discussed, this includes not only Goldstone meson physics per se, but also chiral corrections 
to baryon masses [10], and for the same reason, to mesons containing heavy quarks. 

Recently, also attempts have been made to compute pion and nucleon scattering lengths 
[23,24] from quenched lattice QCD. If one tries to calculate the 1 = pion scattering 
amplitude in quenched ChPT, one actually finds that the imaginary part is divergent at 
threshold, even for nonvanishing pion mass [25]! Again, this can be related to double pole 
terms in the rf propagator. An additional reason is that apparently euclidean quenched 
correlation functions in general cannot be analytically continued to Minkowski space-time. 
(The euclidean four pion correlation functions are well defined.) Further work is needed 
on pion scattering lengths. 
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